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Abstract 

The old idea that an infinite dimensional dynamical system may have its high modes or 
frequencies slaved to low modes or frequencies is re-visited in the context of the 3 D Navier- 
Stokes equations. A set of dimensionless frequencies {ft m (t)} are used which are based on 
L 2m -norms of the vorticity. To avoid using derivatives a closure is assumed that suggests that 
the (l rn (m > 1) are slaved to (the global enstrophy) in the form (l rn = This 

is shaped by the constraint of two Holder inequalities and a time average from which emerges a 
form for T rn which has been observed in previous numerical Navier-Stokes and MHD simulations. 
When written as a phase plane in a scaled form, this relation is parametrized by a set of functions 

1 < A m (r) < 4, where curves of constant A m form the boundaries between tongue-shaped 
regions. In regions where 2.5 < A m < 4 and 1 < A m < 2 the Navier-Stokes equations are shown 
to be regular: numerical simulations appear to lie in the latter region. Only in the central region 

2 < A m < 2.5 has no proof of regularity been found. 

Dedicated to the memory of David Broomhead (1950 - 2014) 


1 Introduction 

1.1 Historical background 

A generation ago a recurrent theme in studies in infinite dimensional dynamical systems was the idea 
that a small subset of low modes or coherent states might conceivably control the dynamics by slaving 
the higher modes to this subset. Having originally emerged from earlier work on centre manifolds in 
studies in ordinary differential equations (Broomhead, Indik, Newell and Rand (1991), Guckenheimer 
and Holmes (1997), Holmes, Lumley and Berkooz (1996)), such ideas obviously have a lasting appeal, 
particularly for those who work in turbulent flows where the number of degrees of freedom are so 
large that resolved computations at realistic Reynolds numbers are hard to achieve: see Moin and 
Mahesh (1998), Donzis, Young and Sreenivasan (2008), Pandit, Perlekar and Ray (2009), Ishihara, 
Gotoh and Kaneda (2009), Kerr (2012, 2013) and Schumacher, Scheelb, Krasnov, Donzis, Yakhot 
and Sreenivasan (2014). For partial differential equations there also emerged a parallel and closely 
related body of work on global attractors and inertial manifolds which aimed to prove the finite 
dimensionality of the system in question in some specified sense: Foias, Sell and Temam (1988), 
Titi (1990), Foias and Titi (1991), Robinson (1996), Foias, Manley, Rosa and Temam (2001). Some 
success was achieved for the one-dimensional the Kuramoto-Sivashinsky equation where an inertial 
manifold was proved to exist (Foias, Jolly, Kevrekidis, Sell and Titi 1988). Further success was also 
achieved for the 2D incompressible Navier-Stokes equations when a global attractor was shown to 
exist with a sharp estimate for its dimension (Constantin, Foias and Temam 1988), with further 
estimates on the number of determining modes and nodes: see Foias and Prodi (1967), Foias and 
Temam (1984), Jones and Titi (1993), Olson and Titi (2003) and Farhat, Jolly and Titi (2014). 

The aim of this paper is to revisit the low/high mode idea in a new way in the context of the 
open question of the global regularity of the 3 D Navier-Stokes equations 

dtu -F u ■ Vm = vAu — Vp + /(*), divw = div / = 0 , (1.1) 

with periodic boundary conditions on a cube of volume V = [0, L]^ er . The body force f{x) is 
taken to be L 2 -bounded and centred around a forcing length scale l in the manner described by 
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Doering and Foias (2002) : in this paper l is taken as t = L for convenience. The two dimensionless 
numbers corresponding to the forcing and the system response are respectively given by the Grashof 
and Reynolds numbers 


Gr = 


i 3/2 ii/ih 


Re = 


LU n 


( 1 . 2 ) 


where Uq = L 3 (Utility with || • ||2 representing the L 2 -norm and (-) T a time average up to time 
T > 0. 


1.2 Some recent analytical and numerical scaling results 

Given the open nature of the question of global regularity of solutions of the 3 D Navier-Stokes 
equations, analyses have conventionally been based on an assumption of some type. The main one 
has been that the velocity field is assumed to remain bounded in some function space, the best 
being the u € L 3 (V) result of Escauriaza, Seregin and Sverak (2003). An alternative computational 
approach has been to discuss low modes in terms of coherent states by projecting onto special 
selections of Fourier-Galerkin modes (Broomhead, Indik, Newell and Rand (1991), Holmes, Lumley 
and Berkooz (1996)). 

Here we break with both of these traditions and instead introduce a set of time dependent 
frequencies or inverse time scales (or ‘modes’) 

(fii(f), n 2 (t), ... n m (t)} , (i.3) 

based on L 2m -norrns of the three-dimensional vorticity field x , t ) which obeys 

(dt + u ■ V) uj = vAld + uj ■ X7u + curl f . (1-4) 

The Q m in (11.31) are defined such that each has the dimension of a frequency 

/ r \ V 2m 

n m (t)= Il~ 3 j^\u\ 2m dv) . (i.5) 

lli(t) is the global enstrophy and the f l m (t) are higher moments. The set (|1.3I) would be widely 
spread if the vector field t ) is strongly intermittent, whereas they would be squeezed closely 
together if u is mild in behaviour. Multiplication by the inverse of the constant frequency wq = vL~ 2 
produces a non-dimensional set = WQ 1 Q rrL . 

In earlier work, a scaled set of the 

~ 2 

An = with a m = - , (1-6) 

4777 — 6 

was shown to have bounded time averages for 1 < m < oo (Gibbon 2011) 

(D m } T < cRe 3 + O (T” 1 ) . (1.7) 

Given that ai = 2, the first of these, (Di) t < cRe 3 , is just the well-known result that the time 
average of the global enstrophy is a bounded quantity. The origin of the a m -scaling in (11.61) comes 
from symmetry considerations. 

It was observed in Donzis et al. (2013) that time plots of the D m (t) from several different 
simulations were ordered on a descending scale. In a further paper (Gibbon et al. 2014), it was 
observed that plots of the maxima in time of In D m /In D\ led to the relation 

D m <Df m ' x with A m \ = (m ~ 1)A +1 , (1.8) 

4 m — 3 
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where the accuracy of the fit lay within 5%. The corresponding fixed values of the fitting parameter 
A lay in the range 1.15 < A < 1.5. These numerical simulations were: (i) a 1024 2 x 2048 decaying 
calculation with anti-parallel initial conditions at about Re\ ~ 400 based on work reported in Kerr 
(2012, 2013); (ii) a forced and a decaying (512) 3 pair of simulations at about Re\ ~ 250 - see Gibbon 
et al. (2014); (iii) data from a large-scale statistically steady simulation (4096) 3 simulation on 10 5 
processors at Re\ ~ 1000, reported in Donzis et al. (2008, 2010) and Yeung, Donzis and Sreenivasan 
(2012). In addition to these, (1 1 .8|f has also been seen in a set of 3H-MHD simulations in similar 
circumstances to those performed for the Navier-Stokes equations: see Gibbon et al. (2015). 


2 A time-dependent closure avoiding derivatives of u> 


The inequality in (11.81) is a numerical observation. To put this, or something close to it, on a rigorous 
footing for the 3D Navier-Stokes equations, the first step is to assume that a strong solutior|j exists 
on a maximal time interval [0, T*) on which it is assumed that fli is bounded - hence all the other 
Cl m are also bounded. Then, instead of taking the standard route of a Sobolev inequality, which 
involves derivatives of u>, the next step is to postulate a relation between Q m and lower Q n which 
could be thought of as a high/low frequency closur^l 

— Tm(^l5 H 2 ? ■ ■ ■ ■ f2 n ) , 1 ^ Tl <C Tfl . (2.1) 


This closure, which is designed to avoid the introduction of derivatives in uj, must be constrained 
and shaped by the fact that Q m must not only satisfy Holder's inequality 

fit < ... < n m < n m+1 ( 2 . 2 ) 


but it must also satisfy a triangular version Holder's inequality for m > 1 (see Appendix A) 

tf) £ (tt) ■ (2 - 3) 

It is well known that the existence and uniqueness of solutions depends entirely on the bounded of 
the 7/i-norm of the velocity field (Leray 1934), which is proportional to . Thus we simplify the 
dependency of Q, m to the first frequency Qi such that our 'closure' in (12.ip is reduced to 

n m = n 1 F m (n 1 ). (2.4) 

The inequalities ()2.2|) and (I2.3[) demand that T m must satisfy both 

Fm+l > JT 2 /^ 2 - 1 ) , Fm> 1. (2.5) 


The first inequality in (I2.5j) can be further simplified by the substitution 


T = 

J m 




m — 1 
m 


( 2 . 6 ) 


for a sequence of smooth, arbitrary functions h m (^fli,rj which must then form a monotonically 
increasing sequence 

1 < h m ^i,rj < h m+ i . (2.7) 

J We are using the standard contradiction method in PDE-analysis where it is assumed that there exists a maximal 
interval time [0, T*) on which solutions of the Navier-Stokes equations exist and are unique: the strategy is to then 
attempt to prove a contradiction in the limit t —► T*. 

'^Strictly speaking, (12.411 is not a closure in the conventional sense used in turbulence modelling. Rather, it is an 
expression of how can be related to lower 12 n , under the constraints (12.2p and (12.311 . without resorting to derivatives 
of uj. Nevertheless, we will continue to use the word ‘closure’ for convenience. 
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In (12.61) and ()2.7p . r is a dimensionless time t = Finally, (12,4j) is reduced to 


^ m — ^1 


hm ( F2 1; T 


( 2 . 8 ) 


The result depends only on the finiteness of the domain and the inverse box frequency zn 0 = uL~ 2 to 
create the dimensionless time r. The relation (12.8p is no more than an expression of the potentially 
arbitrarily large distance between Cl m and in the form of an equality. In fact, ()2.8p contains no 
3 D Navier-Stokes information but once it is considered in this context, the finiteness of the time 
averages in (11.71) must be enforced. An application of a Holder inequality to (12.81) shows that the 
h m must therefore be constrained by 

(C /3 ) t <oo. (2.9) 

One choice of h m consistent with this is that it cannot be any stronger than a power law in Oi 

= l + c m n\ m{T) ~ l ( 2 . 10 ) 

where the monotonically increasing^ A m must lie in the range 

l<A m (r)<4. (2.11) 


Morover, the dimensionless constants c m will be chosen more specifically later. The final result is 
that if there exists a strong solution on the interval [0, T*), then with the choice of h m as in (12.101) . 
there exists a set of exponents {A m (r)} such that 


“ - U) — - - I 




A m (r) —1 


1 < m < OO . 


( 2 . 12 ) 


In practice Oi is so large that we may safely ignore the factor of unity in fi m (Oi,r) in (12.121) . With 
the definition 




1 + (A m (r) — 1) 


771 — 1 


m 


(m - l)A m (r) + 1 

4 m — 3 


(2.13) 


we have 

D m = . (2.14) 

which is a version of (1 1. 8 p with \ m = A m (r) as a function of time. The introduction of the set 
{A m (r)} means the dynamics are analyzed in terms of D\ and {A m (r)} instead of D\ and (L> m (r)}. 
In (1 1. 8 p the inequality is appropriate because in Gibbon et al. (2014) A m was estimated as a constant 
parameter corresponding to maxima in time in plots of \nD m /\nDi versus time. Note that when 
A m = 4, A m ^4 = 1, at which point the D m versus D\ relation is linear. Beyond the range in (12.lip , 
when A m > 4, (12. 141) is no longer valid. Clearly, a further set of numerical experiments are needed 
to analyze more closely the trajectories of A m (r) in (12.141) . 

4 ln Gibbon et al. (2014), in which only the maxima in time of D m were considered, there were indications that the 
Am decreased with m, albeit very weakly, which is not consistent with the constraint Am < A m +i that stems from 
(12.51) . This suggests that the choice of h m made in 112.101) needs a slight modification. This issue needs testing with a 
proper set of numerical experiments. 
































NS-lam4.tex 5th/ll/2015 


5 



Figure ll Plots of the convex curves versus for A m = 2, 3 lying above the straight line A m = 1. As A m increases the 
Clm spread more, corresponding to greater intermittency in the cj-field. 


Am — 4 



Figure 2: A cartoon in the D\ — Dm plane parametrized by Am : in the region 2.5 < A m > 4 the Dm undergo no more 
than algebraic growth with a restriction on (large) initial data for 2.5 < A m < 4; in the region 2 < A m < 2.5 only Leray’s weak 
solutions are known to exist with the only control over Dm being the bounded long time averages ( D m )t < oo; at A m = 2 
the bound on D\ grows exponentially in time; in the region 1 < Am < 2 the Dm lie within an absorbing ball. The dotted 
curves represent the approximate regions of the maxima in the numerical simulations reported in Donzis et al. (2013) and Gibbon 
et al. (2014). The region below the curve A m = 1 is forbidden as Holder’s inequality is violated there. 


3 Navier-Stokes results in the D m — D\ plane parametrized by A m (r) 

The essential idea behind (12.1311 has been to avoid the use of derivatives of u and instead transfer 
the dynamical relationship behind (l rn and fix to the {A m (r)}. This allows us to treat the D m — D\ 
phase as a phase plane parametrized by A m (r), as in Fig. 2. The curves of constant A m are 
drawn as concave curves, although actual orbits A m (r) could potentially wander across the phase 
plane and over these labelled boundaries, even though numerical experiments so far have found that 
their maxima lie in the lower half of the lowest tongue labelled by the dashed curves. We wish 
to demonstrate certain results about the nature of solutions in different sectors of this phase plane 
bounded by curves A m = const, so in this section A m is treated as a constant. 

The curve A m = 1 is associated with the lower bound fli < Cl m , which translates to 

£>“ m/2 < D m where A m l = ot m /2 = m . (3.1) 

4 m — 3 

Thus the regions below the line in Fig. 1 and below the lowest curve in Fig. 2, both corresponding 
to A m = 1, are forbidden by Holder's inequality. Moreover, the relation between D m and D\ is 
linear when X m = 4. The results of i|3.1l alone can be found in Gibbon et al. (2014) : the rest of the 
material is new. 
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3.1 The region 1 < A m < 2 and the curve X m = 2 

From the definition of the D m in (11.61) note that D\ = Lv~ 2 ||u;|| 2 . A purely formal differential 
inequality for D\ is 


\Di<Lv~ 2 1^-v J |Vw| 2 dV + J \Vu\\u>\ 2 dV + L-'WwWzWfWz 

Dealing with the negative term first, an integration by parts gives 

1/2 

( / \u\ 2 dv) 


M 2 dV < 


tv 


\ 1/2 

\Vu\ 2 dV ) 


'V 


>v 


where the dimensionless energy E is defined as 


E = zW 2 //” 1 


I u\ 2 dV. 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


This is always bounded such that 

limt-tooE < cGr 2 . 

Next, the nonlinear term in (13.21) needs to be estimated. The standard result using a Sobolev 
inequality produces a cubic nonlinearity D\ that is too strong for the negative term : all that can be 
deduced from this is that D\ is bounded from above only for short times or for small initial data. The 
difficulty caused by this term has been known for many decades: see Constantin and Foias (1988) 
and Foias et al. (2001). We circumvent this problem by proceeding as follows: 

/ O f 2m—3 1 

/ \Vu\\ui\ z dV = / |a;| m ~ 1 |u>| m ~ 1 \Vu\dV 


< 


< G n 


M 2 dV 


2m—3 

2(m—1) 


<jj\ 2m dV 


2 m (m — 1) 


\Vu\ 2m dV 


1 

2m 


\u\ 2 dV 


'V 


2m —3 
2(m— 1) 


(3.6) 


\u)\ 2rn dV 


2(m— 1) 


2m—3 4m—3 

C T 3 __3 j~\ 2m—2 7~) 2m—2 


1 < m < oo . 


The penultimate line is based on ||Vti|| p < c p ||cd|| p , for 1 < p < oo. Inserting the depletion 

Dm = ci )m T> 1 4m ’ Am gives 

Lv~ 2 [ \Vu\\co\ 2 dV <C2, m tJ7 0 £>i m ’, (3-7) 

Jv 

where £ m ' s defined as 

A m ,\m ( 4m - 3) + 2m - 3 


s m — r,/ \ 

2{m — 1) 

It can now be seen that by using (11.81) . the m-dependency cancels leaving^ 


(3.8) 




1 + ±A, 


(3.9) 


m does not reach its conventional value of 3 unless A m = 
becomes (r = w^t) 


i dD 1 
5 dr 


< -- 


£± 

E 


+ C, m D 


l+^A m 

1 


4. The differential inequality (13.21) now 

+ GrD\ / 2 . (3.10) 


’Lu and Doering (2008) showed numerically that by maximizing the enstrophy subject to divu = 0, two branches 
of the nonlinear term appear, the lower being D i' 78 and the upper Df'" 1 . Later, Schumacher, Eckhardt and Doering 
(2010), suggested that 7/4 and 3 were the likely values of these two exponents: the exponent 1 + |A m — 7/4 
corresponds to A m = 1.5 which lies at the upper end of the range 1.15 < A m < 1.5 observed in Gibbon et al. (2014). 
















NS-lam4.tex 5th/ll/2015 


7 


Given that E is bounded above, D\ is always under control provided X m is restricted to the range 
1 < A m < 2. Formally, there exists an absorbing ball for D\ of radius 

lim^oo-Di < C2 } mGr 2 ~ x m + O (^Gr 4/3 ^j , (3.11) 

It has been shown in Gibbon et al. (2014) that this gives rise to a global attractor A whose Lyapunov 
dimension has been estimated as 


d L (A) < 


C 4 ,m-Re 5 V 2 - 


3 7 6—Ay, 


c 5, 




3 ( 4-Ay, 


(3.12) 


depending on whether one chooses to use the Reynolds or Grashof number. In the limit A m —f 2 the 
radius of the ball in ()3.1ip grows but, at A m = 2, the finiteness of f* Di(t)cIt means that Di(t) 
has an exponentially growing upper bound. 


3.2 The regions 2.5 < A m < 4 

It has also been shown in Gibbon (2012) and Gibbon et al. (2014) that the D m satisfy the differential 
inequality for 1 < m < oo 


Dm < I-K 


+ ro 2,m^ + VO^mGrDln 1/oin 


(3.13) 


where r] m = 2m/3(m — 1) and where w\ )m < zu 2 . m are constant frequencies. The last (forcing) 
term, is hard to handle in conjunction with the others so this will be separated and dealt with last: 
no more than algebraic growth in time can come from it. Let us now divide (I3.13jl by D 3 Z to obtain 

— w l, r n^ rn {t')D m — VJ2 ,m (3-14) 

where 

= ' (3 ' 15) 

A bound away from zero of the time integral of X m (t) is required to show that H“ 2 (f) never passes 
through zero for some range of initial conditions. To achieve this we introduce the nonlinear depletion 
as in (|1.6|) . Noting that r] m + 2 = 2(4m — 3)/3(m — 1) = 2ry rn a“ 1 , it is found that (c m = c^™) 


-y _ ~ T^\Vm (^m 1)(?7T 1)/ TTL 

^-m — G r n U | 

_ ~ i-)2(A m —1)/3 

— 


(3.16) 


It is at this point that we introduce the lower bounc0 on D\(t r ) dt' found 

( 2002 ) 



dt' >tGr + 0(t x ). 


in Doering and Foias 


(3.17) 


°Doering and Foias (2002) have shown that there are two estimates for the lower bound to the integral in (13.171) . 
The first is proportional to tGr and the second to tGr 2 Re~ 2 depending on the relative sizes of the forcing length 
scale t and the Taylor micro-scale. The first has been derived using a Poincare inequality so it is likely to be less sharp 
although the two coincide when the bound Gr < cRe 2 is saturated. For simplicity we use the Gr-bound. 
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This comes into play only if 2(A m — l)/3 > 1, in which case A m > 5/2, where we can then use a 
Schwarz inequality to obtain 

2(A m -l)/3 


rt / rt \ z \ A m J -)/3 

J X m (t') dt! > Cm (J D\ dtj t l-2(A m -l)/3 


> t c m Gr 2(Am_1)/3 + O 


(13.141) can be solved to give 

\[D m {t)] 2 < 


ex 


p{-roi ,rn Jo X rn (t r ) dt'} 


[Dm( o)] 2 - zu 2 , m f 0 exp{—zu lttn f 0 X m (t") dt"}dt' 

exp{-c m tui im tGr 2 ( Am 1 )/ 3 } 

1 
2 


(3.18) 


[D m ( 0)] 2 - w 2 ,m [cm^i, m Gr 2 ^ 4 )/ 3 ] 1 (l - exp{-c m zz7i im t Gr 2 ( A ™ 4 )/ 3 }) 


(3.19) 


(13.191) cannot develop a zero in the denominator if 

1 /2 

D m (0) < Gr (Am_1)/3 , (3.20) 

for any 5/2 < A m < 4, in which case the solution decays exponentially. This initial data is not small 
but it is not huge either. The estimates above have been achieved by neglecting the forcing term in 
(1 3.13 j) . The effect of this term on its own yields 


D m (t ) < [a m ro 3jm Gr (t 0 + t)} a ‘ 


(3.21) 


3.3 The region 2 < A m < 2.5 

The dynamics in the middle tongue-shaped region bounded by the two central curves X m = 2 
and A m = 2.5 in Fig. 2 remains open. Neither the depletion of nonlinearity nor the increase in 
dissipation afforded by (12.141) are strong enough so we must fall back on the existence of Leray's 
weak solutions.The lower bound of Doering and Foias (2002) on the time integral of the enstrophy 
can only be used on the time integral (see (13.181) 

D 2(\ m -i)/3 dt , ( 3 22 ) 

when A m > 2.5. For the range 2 < A m < 2.5 the dimensionless energy E can instead be used to 
bound (13.221) below, but if it passes close to zero for long enough, in the manner of a homoclinic 
orbit, then the resulting lower bound may be too small to be of use. However, ()3.22p could be 
estimated numerically to monitor its behaviour. 



3.4 The region A m > 4 

The relation between D rn and D\ expressed in (12.141) is valid only for 1 < A m < 4. At X rn = 4 the 
relation is linear; namely D m = c m Di. For the region 

D m ^ c m Di ( 3 . 23 ) 

we can insert this relation directly into (13.131) and, provided the c m in (12.141) are chosen such that 
c m > ro 2im ro^, (I3.13P reduces to 

Dm < zuo. m GrD' m l/am , ( 3 . 24 ) 

yielding the same bound for D m (t ) as in (13.21 j) . 
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4 Conclusion 

Subject to the closure relation (j 2.14p . in which the effect of the higher scaled norms D m is hidden in 
the evolution of the set of functions A m (r), it has been shown that the 3 D Navier-Stokes equations 
are regular in each of the regions when X m > 1, with the exception of the region 2 < X m < 2.5. 
There is also the proviso that initial data is limited in the range 2.5 < A m < 4. For initial data set 
in the region X m > 4 (above the straightline in Fig 2), the dynamics are no worse than algebraic 
growth, although such initial data is highly pathological. These results are summarized in Fig. 2. 
The three regular regions are fundamentally different. Solutions lying in the region X m > 2.5 seem 
moribund in the sense that the forcing dominates only algebraically. However, solutions in the lower 
region or tongue 1 < \ m < 2 live in an absorbing ball, the radius of which is given in (13. llh . and 
it is here where all the interesting dynamics lies. It has been shown in Gibbon et aI (2014) that 
solutions here have a corresponding spectrum that is consistent with statistical turbulence theories: 
see Frisch (1995) and Doering and Gibbon (2002). As drawn in Fig. 2 (dotted curves), the large- 
scale numerical simulations reported in Gibbon et al. (2014) lie well within this region. This is only 
partially satisfactory (see below) in the sense that the existence of an absorbing ball is enough for the 
existence of a global attractor provided the solution trajectory A m (r) remains in this region. There 
are two alternatives: 

1. Orbits that originate in the range 1 < A m < 2 always remain there; 

2. Orbits originating in the range 1 < A m < 2 could travel out of this region and into the range 

2 < X m < 2.5 and beyond. However, the nature of this transition is uncertain, and it is unclear 
what the nature of weak solutions would mean numerically if this happened. 

The numerical simulations performed so far have all had their initial data resting in 1 < A m < 2 and 
have shown no evidence for the behaviour in item 2. In fact, the observed range 1.15 < X m < 1.5 
indicates relatively mild dynamics. Unless a rigorous proof is found for the behaviour in item 1, the 
possibility that the behaviour in item 2 could occur for higher values of Re ought to be kept in mind. 

A series of numerical experiments are needed with initial conditions set in the four different regions 

which track the evolution of A m (r) although in the scaled time r = wq t, wq could be so small that 
one may have to compute for significantly large values of the real time t. If the behaviour in item 
2 is observed this would open the question of the physical manifestation of weak solutions. Given 
that these solutions lack uniqueness, would there be there a corresponding physical effect, such as 
multiple branching of the A m -trajectories? 


Acknowledgment 
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A The triangular Holder inequality (I2.3j) 

Consider the definition of 


L 3 Q 2 ™ = f \u\ 2m dV = f \ur\ 2a \<jj\ 2 P dV 

Jv Jv 


>v Jv 

where a + j3 = m. Then, for m > 1 and 1 < p < m — 1 and q > 0, we have 


L 3 n 2 ™ < ( / \u\ ^-p'idv 


\(jj 


1 2 (m+q) c [y 


(A.l) 


oc 

m—p 


m+q 


(A.2) 
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where ^ = i. Solving for a, /3 gives 



q(m-p) p(m + q ) 

a =- and p =-, 

(A.3) 


p + q p + q 

thereby giving 

r>m(p+q) n,q(m-p) 0 p(m+q) 

“m — J c ‘m—p “m +9 

(A.4) 

Now choose q = 

1 and p = m — 1 to obtain 



nZ 2 < , 

(A.5) 


which leads to (12.31) . 
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